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Abstract 

Total reaction cross sections are predicted for nucleons scattering from various nuclei. Projectile 
energies to 300 MeV are considered. So also are mass variations of those cross sections at selected 
energies. All predictions have been obtained from coordinate space optical potentials formed 
by full folding effective two-nucleon (NN) interactions with one body density matrix elements 
(OBDME) of the nuclear ground states. Good comparisons with data result when effective NN 
interactions defined by medium modification of free NN t matrices are used. Coupled with analyses 
of differential cross sections, these results are sensitive to details of the model ground states used 
to describe nuclei. 
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I. INTRODUCTION 



Reaction cross sections from the scattering of nucleons by nuclei (stable and radioactive) 
are required in a number of fields of study; some being of quite current interest [0. An 
example is the transmutation of long lived radioactive waste into shorter lived products using 
accelerator driven systems (ADS). These systems are being designed in the US, Europe, and 
Japan with the added objective of providing an intense neutron source to a sub-critical 
reactor as a new means of energy production and for which nucleon-nucleus reaction cross 
sections are required as well. The technology takes advantage of spallation reactions |2| 
within a thick high-Z target (such as Pb or Bi), where an intermediate energy proton 
beam induces nuclear reactions. The secondary nuclear products, particularly lower energy 
neutrons and protons [^], in turn induce further nuclear reactions in a cascade process. The 
total reaction cross sections of nucleon-nucleus scattering plays a particularly important 
role since the secondary particle production cross sections are directly proportional to them. 
Also they are inputs to intra-nuclear cascade simulations that guide ADS design. As well, 
nucleon-nucleus (NA) cross section values at energies to 300 MeV or more are needed not 
only to specify important quantities of relevance to proton and neutron radiation therapy 
0], but also as they are key information in assessing radiation protection for patients. 

In basic science, these total reaction cross sections are important ingredients to a number 
of problems in astrophysics, such as nucleosynthesis in the early universe and for aspects of 
stellar evolution, especially as the density distribution of neutrons in nuclei are far less well 
known than that of protons. Recently a link has been made between the neutron distribution 
in heavy nuclei, such as 208 Pb, and properties of neutron stars 0, |6| so heightening the need 
for analyses to ascertain the optimal neutron distributions in such nuclei. Analyses of proton 
elastic scattering angular distributions particularly for 200 MeV protons seem to be one such 
method 0. Considering the integral observables of both proton and neutron scattering from 
a given nucleus is another and certainly it will give direct information on the neutron rms 
radius; a property sought in new parity- violating electron scattering experiments |]5], |7|]. 

However, most NA reaction cross sections cannot be, have not been, or are unlikely to 
be, measured. Thus a reliable method for their prediction is required. The usual vehicle 
for specifying these NA total reaction cross sections has been the NA optical potential; a 
potential most commonly taken as a local parametrized function, usually of Woods-Saxon 
(WS) type. However, it has long been known that the optical potential must be nonlocal 
and markedly so, although it has been assumed also that such nonlocality can be accounted 
by the energy dependence of the customary (phenomenological) models |§ . Of more concern 
is that the phenomenological approach is not truly predictive. The parameter values chosen, 
while they may be set from a global survey of data analyses, are subject to considerable 
uncertainties and ambiguities. This is especially true for the optical potentials for nucleon 
scattering from Op-shell nuclei where no such global approach is valid. 

We consider a predictive theory of NA scattering to be one that is "direct" in that all 
quantities required are defined a priori. Thus each result must come from just one run of 
relevant codes, and there should be no post evaluation adjustments save for allowance of 
known a priori uncertainties in the specification of the input information. 

With the nucleus viewed as a system of A nucleons, NA scattering is determined to first 
order by an optical potential formed by folding, with a suitable specification of the ground 
state density of the target, appropriate interactions of that projectile with each and every 
nucleon within the nucleus. Over the past decade or more, such microscopic approaches 
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defining the NA optical potential have been quite successful in predicting elastic scattering 
data (differential cross sections in particular). Studies forming the optical potential in 
both momentum and coordinate space have been made with success ||. In the coordinate 
space approach, and for analyses that are based upon the DWBA programs of Raynal 
H 0IJ, the projectile-target nucleon interaction takes the form of a complex, energy and 
density dependent, effective NN interaction. Of those programs, DWBA98 has been used 
to evaluate all of the cross sections shown later herein. Appropriate effective interactions 
can and have been defined that, upon folding with good structure wave functions of nuclei, 
give credible optical potentials. Using those optical potentials, differential cross sections 
and spin observables such as the analyzing powers for proton scattering at many energies in 
the range 40 to 800 MeV (65 and 200 MeV in particular) and from diverse targets ranging 
from 3 He to 238 U have been predicted and found to have excellent agreement with data. 
Moreover, and very recently @ , it has been shown that analyses of differential cross sections 
of proton elastic scattering can select between alternative model predictions of the neutron 
rms radius in 208 Pb; such are sensitive also to the surface distribution of its neutron matter. 

The microscopically formed optical potentials are complex and energy dependent from the 
like properties of the effective NN interaction ||. Such properties arise from mapping the 
effective interactions to iViV g matrices that are solutions of the Bruckner-Bethe-Goldstone 
(BBG) equations for nuclear matter. The BBG equations carry medium modification of 
the NN scattering due to Pauli blocking and to a background mean field. Details of the 
effective interactions, of the folding process that gives the (nonlocal) optical potential, and of 
the successful predictions found therefrom of differential cross sections and analyzing powers 
from the scattering of protons at diverse energies and from diverse targets, are given in the 
literature [g]. 

II. ELEMENTS OF THE OPTICAL POTENTIAL FOR NUCLEON SCATTERING 

Formally, the nonlocal optical potentials from a first order folding model can be written 

U(ri,r 2 ;E) = jtf(ri - r 2 ) J (f* n (s) v D [r u ) <p n (s) ds + <p*(ri) v Ex (r 12 ) ^„(r 2 ) j 

U D { ri ; E)5(n - r 2 ) + U Ex {r u r 2 ; E) , (1) 

where v E , v Ex are combinations of the components of the effective NN interactions, ( n are 
ground state nucleon shell occupancies (more generally they are the ground state OBDME), 
and <f n {x) are nucleon bound state wave functions; denoted SP functions hereafter. All 
details and the prescription of solution of the associated nonlocal Schrodinger equations are 
given in the review ||. 

The results to be discussed have been found by solving the actual nonlocal Schrodinger 
equations defined with potentials as given (formally) by Eq. ([I]). Two formulations of those 
optical potentials have been used. They and the results are identified by the appellations, 
g- and t- folding, according that the effective NN interactions have been defined by their 
mapping to the BBG g matrices or to the basic free NN scattering t matrices respectively. 
The latter are solutions of Lippmann-Schwinger (LS) equations. In both cases, the driving 
NN interaction has been the Bonn-B NN potential [pTTf] . 

From practical necessity the model descriptions of nuclei in the mass range 3 to 238 vary 
in complexity. With the light mass nuclei (A < 12 for example), quite large and complete 
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shell model spaces with potentials either fitted or formed as G-matrix elements have been 



made JT2, [13], [14|]. While large space model studies of heavier nuclei are being sought, the 
dimensions of the problem preclude our use of all but Ohu shell model specifications for 
most heavier nuclei. Indeed, for targets heavier than mass 90 we have used an even simpler, 
packed shell, definition of their ground states. Nevertheless with such model prescriptions 
and using harmonic oscillator (HO) SP functions with oscillator energies selected according 
to an A~3 rule, very good predictions of the scattering of 65 and 200 MeV protons have been 
obtained for all but the light mass nuclei ||. However, 208 Pb and 40 Ca are special cases. 
Recently |5j, [/J, a Skyrme-Hartree-Fock (SHF) model of those nuclei was made and OBDME 
required in our folding procedure were extracted. The associated density distributions vary 
noticeably from that given by the HO (packed) shell model and, not surprisingly, so do 
proton differential cross sections. 

For nuclei with A < 12 typically, better spectroscopy is needed. So also are more realistic 
matter distributions for studies with light mass exotic nuclei, such as of radioactive beam 
scattering from hydrogen. For example, the reaction cross section for 40. 9A MeV 6 He 



scattering from hydrogen [13 varies from 350 mb, found when He has a neutron skin 
as expected with a standard shell model description, to 406 mb when that distribution 
is extended further to be classified as a halo by choosing valence neutron SP functions 
consistent with the single neutron separation energy in 6 He. The measured value is 409 ± 
21 mb IP, HI. 



III. PHASE SHIFTS, S-MATRICES, AND OBSERVABLES 

Irrespective of the means used to define NA optical potentials, the objective is to de- 
fine the S matrix, or equivalently the (complex) phase shifts Sf(k), where the superscripts 
identify the values j = I ±1/2. These relate by 

Sf{k) = e 2lS ^ = r^(fc)e 20f K ±(fc) ] (2) 

where 

V ±(k) = \S±(k)\=e-^^]. (3) 

With E oc k 2 , the elastic, reaction (absorption), and total cross sections respectively then 
are given by 

__ 00 

a *w = ^J2{( l+1 )\ s ^ k )- 1 \ 2 + l \ s r(k)-i\ 2 } , 



1=0 

oo 

1=0 



and 



a TOT (E) = a c i(E) + a R (E) 

1=0 
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The scattering amplitudes are then 2x2 matrices in the nucleon spin space and have the 
form, 



f(6)=g{0) + h{6)tT-h. (4) 

where 

9( ) = I E W + X ) ^(*0 ~ !] } P ^ ■ 

z=o 

= i£[#(*)-^(*)l (5) 



2=1 

In terms of these (complex) amplitudes, the (elastic scattering) differential cross section is 
defined by 

d° — ,2 . , , /„m2 



and the analyzing power A y {9) by 



^)r + W)r . (6) 



29K[g*(fl)fc(0)] m 



IV. RESULTS OF CALCULATIONS 



All results we show have been evaluated using the DWBA98 program |L0| , input to which 
are density dependent and complex effective NN interactions having central, two-nucleon 
tensor, and two-nucleon spin-orbit components. The effective interactions we use have been 
generated for energies from 10 MeV to over 300 MeV in 10 MeV steps by an accurate mapping 
to iViV t- and g matrices found by solutions of the LS and BBG equations respectively and 
based usually upon the Bonn NN potentials. Details are given in the review ||. 

Other input to DWBA98 are the ground state occupancies (or OBDME) and the associ- 
ated SP functions. The SP functions used in most of the calculations for nuclei of mass 20 
and above at best come from a Ohu shell model which has been adopted to describe their 
ground state occupancies. For the lighter mass nuclei considered, larger shell model spaces 
were used to define their ground states, and in some cases the interaction potentials defined 



as G matrix elements of a realistic interaction fT2|. In shell model studies using those G 



matrix elements, the oscillator energy {Tiuj) for the SP functions also is specified. As stated 
earlier, the cases of 208 Pb and 40 Ca are special in that we have used structure information 
taken from recent SHF studies [H, 17]]. 



A. Energy variation of proton total reaction cross sections 

In this subsection we present our predictions of the total reaction cross sections for proton 
scattering up to 300 MeV for diverse nuclei, ranging in mass from 6 Li to 238 U. In all cases, 
at least two calculations were made. The first of these used the effective interaction defined 
from the t matrices of the Bonn-B interaction while with the second, that built upon the 
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TABLE I: Data source table for proton and neutron reaction data used 



Nucleus 



Proton references (in year order) 



6 Li 
9 Be 

12 C 
16Q 
19 F 

27 Al 

40 Ca 
63 Cu 

90 Zr 

118 Sn 

140 Ce 
i59 Xb 

181 Ta 
197 Au 

208p b 
238 TJ 
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associated g matrices was used. Comparison of the results of each pair of calculations 
demonstrates the effects in predictions due to the medium modification to the free NN 
interaction that define the g matrices. The ensuing t- and g-folding results are portrayed in 
the figures by the dashed and solid curves respectively. 

As noted, the structure models of the light mass nuclei involve diverse complete NTiuo 
bases. For the lightest, 6 Li, a complete (0 + 2 + 4)hu model of structure has been used, 
while for 9 Be and 12 C the OBDME have been defined from complete (0 + 2)hu shell model 
calculations |17|. In addition we have calculated the reaction cross sections from 118 Sn and 



159 Tb allowing the outer (neutron) shell to have a smaller (15-20%) harmonic oscillator 
energy. By that means, the neutron surface of each is slightly more extended than with 
the base (packed shell) model forms; a very simple allowance for any effect of ground state 
deformations. This idea for varied surface SP functions within an HO model has been used in 



the guise of a two frequency shell model |18| . Our wave functions are not so well determined 
of course. 

The results for scattering from 6 Li, 9 Be, 12 C and from 16 are displayed respectively in 
segments (a), (b), (c) and (d) of Fig. |l[ The experimental data shown therein were taken 
from the references listed against those nuclei in the data table, Table |; which also lists 
the sources of the data that are shown in the figures to follow. The data from the four 
lightest mass nuclei are well reproduced by g-folding calculations made with large space 
shell model structure (solid curve), but they are not with t-folding calculations (dashed 
curve). However, large space structure calculations are necessary if one is to describe the 
physics even adequately. As indicated above, for 6 Li, such structure was found from a shell 
model calculation []5D| made using a compete (0 + 2 + 4)huj space while those for 9 Be, 12 C 
and 16 were made using complete (0 + 2)hu models. Results found using the simpler Ohu 
structure model of Cohen and Kurath pT| within g-folding, underestimate the data at most 
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FIG. 1: Energy dependencies of cfr for proton scattering from (a) 6 Li, (b) 9 Be, (c) 12 C, and (d) 
16 0. The solid curves are predictions made using ^-folding optical model calculations while those 
found using i-folding optical potentials are displayed by the dashed curves. 



energies, and for 9 Be particularly. This reflects the too compressed density profile for the 
nuclei given by the simpler model. Results are displayed for proton energies from 10 MeV. 
Although experimental data exist to lower energies in these cases, we do not consider the 
first order folding prescription for the optical potential to be appropriate in the low energy 
regime of scattering from these nuclei. Up to and over 20 MeV excitation, their spectra have 
numerous distinguishable states. Excitations to regions of low level density are not taken 
into account forming the optical potentials. We have confidence in the optical potentials 
when the input energy coincides to excitations to regions of high level density and where 
particle emission is feasible. 

The 12 C results in particular are worth comment. For this nucleus, as with the others, 
the reaction cross sections obtained from those ^-folding calculations for 12 C are in very 
good agreement with the experimental data up to 300 MeV. Most evidently, the medium 
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effects differentiating the g- from the t- matrices used in the folding scheme defining the 
optical potentials are required for predictions to match observation. The t-folding model 
overestimates the data by 20-40% within the energy regime below 200 MeV. But some data, 
notably at 61 MeV |J^] and at 77 MeV MeV with 12 C, are in disagreement with both 
calculated results. 

Predictions for proton scattering from 19 F, 27 A1, 40 Ca, and from 63 Cu respectively are 
compared with the data in the segments (a), (b), (c) and (d) of Fig. ||. (As noted above 
the wave functions for 40 Ca were obtained from an SHF model j7|.) Again all data are well 
reproduced by the ^-folding calculations. With 27 Al however, three data points between 180 
and 300 MeV are in better agreement with the results of t-folding calculations while one 



data point, at 61 MeV |32| , is in disagreement with both calculations. 

Data are shown in this figure again from 10 MeV but with 40 Ca in particular, the folding 
model approach is not expected to be reliable at the energies in the range 10 to 20 MeV. 
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FIG. 3: Energy dependencies of gr for proton scattering from (a) 90 Zr, (b) 118 Sn, (c) 140 Ce, and 
(d) 159 Tb. Basic notation is as for Fig. |l|, but now also with dotted curves that display predictions 
obtained with the g-folding model potentials formed with the outer most shell neutrons specified 
by a harmonic oscillator with an oscillator length 10% larger. 



Those excitations energies correspond to a region of low level density in Ca. Indeed the 
reaction data from 40 Ca show rather sharp resonance- like features below 20 MeV. For 63 Cu 
however, no such sharp structures are evident in the reaction cross section data and our 
prediction with a g-folding potential at 10 MeV gives a value in quite reasonable agreement 
with observation. With both 40 Ca and 63 Cu, the g-folding results are in very good agreement 
with the data for energies above 20 MeV. That is in stark contrast to the t-folding results. 
The t-folding results underestimate the data below 20 MeV and overestimate considerably 
the data above 40 MeV. In the 20 to 40 MeV zone, both calculations give results in reasonable 
agreement. Such trends are evident for most heavy nuclei. 

In Fig. H we present the data and our predictions of the total reaction cross sections for 
proton scattering from 90 Zr, 118 Sn, 140 Ce and 159 Tb. They are shown in segments (a), (b), (c), 
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and (d) respectively and compared against data taken from the relevant references given in 
Table |. Again results from ^-folding calculations are in very good agreement with the data 
while the t-folding results are overestimates at and above 40 MeV and underestimates the 
data below 20 MeV. The third p- 118 Sn total reaction cross section result given in segment (b) 
of Fig. |3] and portrayed in that figure by the dotted curves, was obtained from a g-folding 
optical potential formed by varying the surface neutron orbit (/in/2) to be that for an 
oscillator length increased by 10% from our basic calculation. With the (slightly) extended 
neutron distribution that results, the ^-folding potential total reaction cross sections then 
are in very good agreement with the data; save for the ubiquitous 61 MeV value. Likewise 
there is a datum at 32 MeV at odds with our results. But that point also is at odds with 
other data. Also in the mismatch collection, the datum at 17.5 MeV p- 140 Ce scattering is 
underestimated. For 159 Tb, the (/-folding result (solid curve) is still a quite good replication 
of data but the calculations obtained from g-folding optical potentials formed by varying 
the surface neutron orbit (hg/2) to be that for an oscillator length increased by 10% (dotted 
curve) are better. 

In segments (a), (b), (c) and (d) of Fig. f|, we compare the calculated total reaction cross 
sections with proton scattering data from 181 Ta, 197 Au, 208 Pb, and 238 U. Again all data are 
best and well described by the predictions we make with the ^-folding model. There is a cross 
over regime in energy where the predictions of our t- and ^-folding potential equate. The 
differential cross sections do not however and when such are also considered ||, preference 
is given to the g-folding potential results at all energies. 

The dotted curve in the case of 208 Pb [segment (c)] results when the oscillator length for 
the outer neutron shell (213/2) of the simple packed shell model we have used to describe 
the nucleus is increased by 8%. The associated increase in the matter profile brings the 
predicted reaction cross sections then in very good agreement with observation. Using the 
SHF wave functions 0, gives the result displayed by the dot-dashed curve in this figure. 
Clearly using these new functions has made a slight change to the predictions found with the 
simple HO packed model (solid curve); the modified HO model result is in better agreement 
with that of the SHF model. 

Exceptional data points are found at 19.8 MeV in 181 Ta where that data point is un- 
derestimated with our calculation by 20%, at 29 MeV in 197 Au, and in 208 Pb near 30, 61 
and 77 MeV. Such exceptional points were also noted in the data from lighter mass nuclei. 
But in most cases, those exceptional point values do not agree with other measurements 
made at close values of energy. For example other data from 208 Pb taken at 60.8 MeV |36 



and at 65.5 MeV |28| give different results and in fact reaction cross section values that are 



consistent with our predictions. We note that Menet et al. [36] argue for a much larger 



systematic error in the studies reported in the relevant earlier experiments. 



B. Mass variation of proton total reaction cross sections 

The mass variations of total reaction cross sections for the scattering of 25, 30, 40, 65, 100 
and 175 MeV protons are shown in the different segments (as labeled) of Fig. |5[ From that 
figure, it is evident that the ^-folding results are in quite good agreement with data while 
the t-folding results underestimate most of the 25 MeV data, are in reasonable agreement 
with the 30 and 40 MeV data but overestimate most of the 65 MeV data. Note that the 
extended matter SP states have been used with the 208 Pb and the Sn isotopes calculations. 

The disparities between the t- and (/-folding potential results for the reaction cross sections 
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FIG. 4: Energy dependencies of ctr for proton scattering from (a) 181 Ta, (b) 193 Au, (c) 208 Pb, and 
(d) 238 U. Basic notation is as for Fig. |l[ For 208 Pb, the dotted curve is the result of extending the 
iis/2 neutron orbit by increasing the oscillator length for that shell by 8%. The dot-dashed curve 
is the predictions found using SHF wave functions. 



are more evident at higher energies. In segments (e) and (f) of Fig. |5| we display the mass 
variation of the total reaction cross sections measured []52| at 100 and 175 MeV. Again, 
the ^-folding model predictions are in excellent agreement with the measured values, while 
the t-folding results overestimate observations typically by 150 mb. At 100 MeV proton 
scattering, proton total reaction cross sections from many nuclei in the mass range to 238 U 
have been measured and it is very clear that the g-folding model predictions are in good 
agreement with them. Fewer measurements have been made at 175 MeV, but they too 
span the mass range to 238 U and the results of those measurements also are in very good 
agreement with the ^-folding optical model predictions. 

We have seen in most of the previous figures that, at some energy point, the total reaction 
cross sections obtained by the t-folding and g-folding calculations are same. In Fig. |6| we 
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FIG. 5: The mass variations of proton reaction cross sections at diverse energies. Scattering at 
energies of 25, 30, 40, 65, 100, and 175 MeV are presented in segments (a), (b), (c), (d), (e), and 
(f), respectively. 



display those crossing energy values as a function of A 1 ^; the energy is in MeV. It clearly 
indicates that there is a line in mass and energy above and below which the ^-folding 
calculation results are smaller and larger respectively than the corresponding i-folding ones. 



C. Energy/mass variations of neutron total cross sections 

Accurate measurements of neutron total reaction cross sections are far more difficult to 
achieve than their proton counterparts. Indeed usually those cross sections are obtained 
by subtracting the elastic from the total scattering cross section. While both of those 
cross sections can be measured with some accuracy, the subtraction of two large numbers 
with attendant uncertainties is subject to numerical problems. We note, however, that 
a new technique from Japan for measuring neutron total reaction cross sections utilizing 
in-beam and out-beam methods similar to those used in proton scattering shows promise. 
Nevertheless, herein, we concentrate on analyses of total scattering cross section data. 
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FIG. 6: The mass-energy equivalence points for t- and ^-folding results for proton total reaction 
cross sections. The line is drawn only to guide the eyes. 



et al. 



The data we have chosen to analyze have been taken from a recent survey by Abfalterer 
d3[]. That survey includes data measured at LANSCE that are supplementary and 

From that recent data compilation, 
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additional to those published earlier by Finlay et al. 
we have selected cases for study having target masses that span the stable mass range and 
which includes as many of the nuclei as possible for which we have analyzed proton scattering 
data. Specifically we calculate neutron total scattering cross sections for 6 Li, 12 C, 19 F, 40 Ca, 
89 Y, 184 W, 197 Au, 208 Pb and 238 U. In some cases the data taken were from natural targets; 
we indicate those as discussion of them occurs in the text. Given that analyses with the 
g-folding optical potential has been found the more appropriate with the proton scattering 
analyses, we only present in the ensuing figures, results obtained from the g-folding optical 
models for neutron scattering. In forming those neutron optical potentials, the structure 
models used were those that defined proton scattering. The effective A A interactions were 
as well. 

Predictions of the total cross sections for neutrons scattered from 6 Li, 12 C, 19 F, and 40 Ca 
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FIG. 7: Total cross sections for neutrons scattered from (a) 6 Li, (b) 12 C, (c) 19 F, and (d) 40 Ca. 



The data are those of Abfalterer et al. [53]; in (b) the data correspond to scattering from natural 
carbon. 



are compared to the data in Fig. 0. The models for 6 Li, 12 C, and 40 Ca to specify the densities 
were those used in the calculations of proton scattering. The model for 19 F was a Ohu shell 
model calculation using the WBT interaction of Warburton and Brown The oscillator 
parameter for the HO SP functions in that case was 1.855 fm, chosen to reproduce the rms 
radius of 19 F (2.9 fm [56]). In all cases except for 40 Ca, there is excellent agreement with 



the data, and the turnover in the cross sections at ~ 20 MeV for natural carbon and 19 F are 
predicted. The low-energy structure in the cross section for 40 Ca is reproduced, although 
the peak at 40 MeV is underpredicted by 10%. Also for 40 Ca the cross section at 300 MeV 
is underpredicted. 

Our predictions of the total cross sections for neutrons scattering from 89 Y, 184 W, 197 Au, 
and 238 U are compared to the data in Fig. [| In the case of 184 W, the data were taken from 
scattering using natural tungsten as the target. The structure of three of these nuclei was 
obtained from a simple packed model as used previously. For 89 Y, however, the shell model 
of Ji and Wildenthal p7| was used to specify the density. The four results exhibit the same 
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FIG. 8: Total cross sections for neutrons scattered from (a) Y, (b) W, (c) Au, and (d) 38 U. 



The data are those of Abfalterer et al. [53]; in (b) the data correspond to scattering from natural 
tungsten. 



features. The predicted cross sections agree well with the experimental results between 70 
and 200 MeV, with the peak cross section values slightly underpredicted. Below 70 MeV the 
cross sections are overpredicted but the structural character exhibited in the data is found. 
To obtain better agreement at energies 40 to 70 MeV at least, improvements in the SP 
functions are required. Improvements such as by using a SHF model seem needed, given the 
results we shall present next for 208 Pb. Above 200 MeV, all results underpredict the data. 
We believe that this indicates that our present effective interaction for those energies at and 
about the pion production threshold need be improved. Possibly more explicit contributions 
from the A resonance in the NN interaction in defining the effective NN force to be used 
in the g-folding giving the optical potentials are needed . 



As stated above in the case of proton scattering, 208 Pb is a special case. To specify 
the density of 208 Pb we have used the SHF calculation of Brown |5|, PI. The results of our 
calculations of the total neutron scattering cross section from 208 Pb using that model are 
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FIG. 9: Total cross section for the scattering of neutrons from 208 Pb. The results of our calculations 
are compared to the data of Finlay et al. |54|] (circles) and of Abfalterer et al. [^] (squares). The 
latter set correspond to scattering from natural Pb. 



given in Fig. |9|. Therein, we compare our prediction with the data of Finlay et al. [Q and of 
Abfalterer et al. [53|. The latter set correspond to scattering from natural Pb. From 60 to 



200 MeV, the agreement with the data is excellent, as it was for the other nuclei also. Below 
60 MeV, the energy dependence is reproduced, although the minimum at 50 MeV is slightly 
overpredicted. Note that when the oscillator model is used, the calculated results do not 
reproduce the data in this energy regime as adequately 0. Above 200 MeV, the predicted 
cross section falls too sharply as it did for the other scattering cases. The variation of the 
structure model from HO to SHF did not influence much the result above 200 MeV. Again we 
think that our effective force may be at fault for this data at these energies. The mismatch 
above 200 MeV is influenced by the mass of the target however, with little problem evident in 
the light mass results as shown in Fig. ^. Given that those light mass nuclei as characterized 
as mostly "surface" , we also believe that the problems with the effective interactions at the 
higher energies relates to the character of the force at central densities. 
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V. CONCLUSIONS 



A microscopic model of the nucleon-nucleus optical potential in coordinate space has 
been used to predict successfully the total reaction cross sections of nucleons from nuclei. 
That optical potential has been formed by folding complex energy- and density-dependent 
effective NN interactions with OBDME of the target obtained primarily from shell models 
of the nuclei. As the approach accounts for the exchange terms in the scattering process, 
the resulting complex and energy dependent optical potential also is non-local. We have 
found that it is crucial to use effective NN interactions which are based upon realistic free 
NN interactions and which allow for modification from that free NN scattering form due to 
nuclear medium effects of Pauli blocking and an average mean field. For optimum results, 
and for the light masses in particular, it is essential also to use the best (nucleon based) model 
specification of nuclear structure available. Marked improvement in results were obtained 
when, for 9 Be and 12 C in this study, complete (0 + 2)%u shell model calculations were used 
to define the OBDME required in the folding processes. For scattering from heavy nuclei at 
around 300 MeV, our results indicate the need for improvement in the effective NN force 
possibly by explicit inclusion of A effects but also of those associated with the interaction 
at central field densities. 
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